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, Abstract 

Q^' We consider the Dirac equation with a magnetic-solenoid field (the superposition 

. of the Aharonov-Bohm solenoid field and a collinear uniform magnetic field) . Using 

\^ I von Neumann's theory of the self-adjoint extensions of symmetric operators, we 

T-H ' construct a one-parameter family and a two-parameter family of self-adjoint Dirac 

■ Hamiltonians in the respective 2-1-1 and 3-1-1 dimensions. Each Hamiltonian 
^ ! is specified by certain asymptotic boundary conditions at the solenoid. We find 

I the spectrum and eigenfunctions for all values of the extension parameters. We 

■ also consider the case of a regularized magnetic-solenoid field (with a finite-radius 
. solenoid field component) and study the dependence of the eigenfunctions on the 
I behavior of the magnetic field inside the solenoid. The zero-radius limit yields 

■ a concrete self-adjoint Hamiltonian for the case of the magnetic-solenoid field. In 
. addition, we consider the spinless particle in the regularized magnetic-solenoid field. 
I By the example of the radial Dirac Hamiltonian with the magnetic-solenoid field, 

Q^' we present an alternative, more simple and efficient, method for constructing self- 

D . adjoint extensions applicable to a wide class of singular differential operators. 



' 1 Introduction 



In the present Chapter, we study the Dirac particle in a magnetic-solenoid field. This 
field is the superposition of the Aharonov-Bohm (AB) field (the field of an infinitely long 
and infinitesimally thin solenoid) and a collinear uniform magnetic field. We note that 
the solutions of the Schrodinger, Klein-Gordon, and Dirac equations in such a field were 
already discussed in [1, 2, 3, 4]. In particular, the solutions of the Dirac equation with the 
magnetic-solenoid field in 2 -|- 1 and 3-1-1 dimensions were studied in detail in [4] . These 
solutions were used to calculate various characteristics of the particle radiation in such a 
field [5]. In fact, the AB effect in the synchrotron radiation was investigated. However, a 
number of important and interesting aspects related to the rigorous quantum-mechanical 
treatment of the solutions of the Dirac equation with the singular magnetic-solenoid field 
were not considered. In particular, it was pointed out that a critical subspace exists where 
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the problem of self-adjointness for the Dirac Hamiltonian taken naively arises [4]. This 
problem and the associated problem of the completeness of the solutions was not solved. 

We note that even for the pure AB field it was not simple to solve these two problems 
inherent in both the relativistic and nonrelativistic cases. The construction of self-adjoint 
nonrelativistic Hamiltonians for the AB-ficld case by the method of the self-adjoint ex- 
tensions of symmetric operators was first studied in detail in [6] where the regularized 
AB field was also considered. The self-adjoint extension method was also invoked in the 
anyon physics [7, 8, 9]. The self-adjointness problem and the need for self-adjoint extension 
method in the case of the Dirac Hamiltonian with the pure AB field in 2 -|- 1 dimensions 
was recognized in [10, 11]. The interaction between the magnetic momentum of a charged 
particle and the AB field essentially changes the behavior of the wave functions at the 
magnetic string [11, 12, 13]. It was shown that there exists a one-parameter family of 
self-adjoint extensions specified by certain boundary conditions at the origin; in what fol- 
lows, we call such boundary conditions the self- adjoint boundary conditions. Self- adjoint 
extensions for the case of the Dirac Hamiltonian in 3 + 1 dimensions were constructed in 
[14] . An alternative method for solving the Hamiltonian extension problem in 2 + 1 and 
in 3 -|- 1 dimensions was presented in [15, 16]. It was shown that in 2 -|- 1 dimensions, 
only two values of the extension parameter correspond to the presence of the point-like 
magnetic field at the origin, whereas other values of the extension parameter correspond 
to additional contact interactions [17]. One possible self-adjoint boundary condition was 
obtained in [12, 18, 19] by specifically regularizing the Dirac delta function, requiring the 
continuity of the both Dirac spinor components at a finite radius, and then shrinking 
this radius to zero. Other extensions in 2 + 1 and 3 + 1 dimensions were constructed in 
[20, 21, 22] by imposing spectral boundary conditions of the Atiyah-Patodi-Singer type 
[23] (the MIT boundary conditions) at a finite radius, and then taking the zero-radius 
limit. It was shown that for some extension parameters it is possible to find a domain 
where the Hamiltonian is self-adjoint and commutes with the hclicity operator [24, 25]. 
The bound state problem for particles with magnetic moment in the AB potential was 
considered in detail in [26, 27, 28]. The physically motivated boundary conditions for the 
particle scattering by the AB field and a Coulomb center were studied in [29, 30, 31]. 

It is unclear whether the self-adjointness problem for the case of the magnetic-solenoid 
field can be automatically solved by directly extending the results for the pure AB field, 
i.e., by applying the same boundary conditions because the presence of the uniform mag- 
netic field essentially changes the domains of the relevant operators, in particular, the 
deficiency subspaces, and changes the energy spectrum of the spinning particle from con- 
tinuous to discrete. Therefore, the case of the magnetic-solenoid field requires an inde- 
pendent study. By analogy with the pure AB field, it also seems important to consider 
the regularized magnetic-solenoid field (we call the regularized magnetic-solenoid field the 
superposition of a uniform magnetic field and the regularized AB field) and to study the 
solutions of the Dirac equation in such a field. This problem was not solved before and 
is of particular interest irrespective of the extension problem. The Pauli equation with 
the magnetic-solenoid field was recently studied in [32, 33]. The problem of defining the 
Hamiltonian in a particular case of the magnetic-solenoid field (both fields have the same 
direction) was considered in [34] (the scalar case) and in [37] (the spinning case in 2 -|- 1 
dimensions). We note that the AB symmetry is violated for the spinning particle case, 
which is therefore sensible to the solenoid fiux direction. As a consequence some of the 
results can depend on the mutual orientation, parallel or antiparallel, of the solenoid fiux 
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and the uniform magnetic field: either parallel or antiparallel. We study both possibili- 
ties in detail. The 3+1 dimensional spinning problem and the relation of the self-adjoint 
extensions to the regularized problems were not studied. 

In this Chapter we consider the Dirac equation with the general magnetic-solenoid 
field (the uniform magnetic field and the AB field can have both the same and opposite 
directions) and with the regularized magnetic-solenoid field in 2 + 1 and 3 + 1 dimen- 
sions. We start with separating the variables and describing the exact square-integrable 
solutions of the radial Dirac equation in the magnetic-solenoid field in 2 + 1 dimensions. 
Using von Neumann's theory of the self-adjoint extensions of symmetric operators, we 
then construct a one-parameter family of self-adjoint Dirac Hamiltonians specified by 
boundary conditions at the AB solenoid and find the spectrum and eigenfunctions for 
each value of the extension parameter. We reduce the (3 + l)-dimensional problem to the 
(2 + l)-dimensional one by a proper choice of the spin operator, which allows realizing all 
the programme of constructing self- adjoint extensions and finding spectra and eigenfunc- 
tions in the previous terms. We then turn to the regularized case of finite-radius solenoid. 
We study the structure of the corresponding eigenfunctions and their dependence on the 
behavior of the magnetic field inside the solenoid. Considering the zero-radius limit with 
the fixed value of the magnetic flux, we obtain a concrete self-adjoint Hamiltonian corre- 
sponding to a speciflc boundary condition for the case of the magnetic-solenoid field with 
the AB solenoid. For completeness we also study the behavior of the spinless particle in 
the regularized magnetic-solenoid field. 

It is well worth noting that in our particular case, as in many other cases where 
singular differential operators occur, the general von Neumann's procedure for the self- 
adjoint extensions of symmetric operators can be significantly reduced to analyzing the 
behavior of the corresponding functions in the vicinity of the singularity. In Sec. 5, we 
present this more simple method for the self-adjoint extensions as applied to the radial 
Dirac Hamiltonian with the magnetic-solenoid field. 

2 Exact solutions 

We consider the Dirac equation {c — h—l)m3 + l and 2 + 1 dimensions, 

ido^ = H-^, = 7° (7P + M) . (1) 

Here, 7^ = (70,7) , 7 = (7'=), = idk - qA^, A; = 1,2 for 2 + 1 and k = 1,2,3 for 

3 + 1, 1/ = (0, fc) ; q is an algebraic charge, for electrons, q = —e < 0. As an external 
electromagnetic field, we take the magnetic-solenoid field. The magnetic-solenoid field is 
the coUinear superposition of the constant uniform magnetic field B and the Aharonov- 
Bohm field B^^ (the AB field is a field of an infinitely long and infinitesimally thin 
solenoid). The complete Maxwell tensor is 

F,.^B{SlSl-SiSl),B^B^^ + B. 

The AB field is singular at r = 0, 

B^^ = (!>6{x')5{x^) . 

The AB field creates the magnetic fiux $. It is convenient to represent this flux as: 

$ = (io + fi) $0, ^0 = 27r/e , (2) 
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where Iq is integer and < /i < 1. 

If we use the cyhndric coordinates ip and r: = rcosip, = rsin<^, then the 
potentials have the form 

„ . r, . / XT sin(^ , . / M cos (/J 

^0 = 0, eAi ^[lo + ii + A (r)] -y-, eA2 ^ - [lo + + A (r)] , 

(^ = in 3 + 1), A{r)=eBry2. (3) 

We conventionally treat the Dirac equation quantum-mechanically as the evolution 
Schrodinger-type equation with the Hamiltonian H, an operator in the appropriate Hilbert 
space, and restrict ourselves to the stationary solutions, i.e., to the (generalized) eigen- 
functions of the Dirac Hamiltonian. This implies that we seek the stationary solutions 
that are bounded at infinity and locally square integrable. 



2.1 Solutions in 2+1 dimensions 

We first consider the problem in 2 + 1 dimensions There are two non-equivalent represen- 
tations for the 7-matrices in 2 + 1 dimensions: 

/ = a\ 7' = 7' = -^^'C, C = ±1 , 

where the "polarizations" ( = ±1 correspond to the respective "spin up" and "spin down" 
particles and cr = (a*) are the Pauli matrices. The stationary case is assigned the following 
form of the spinors \l/(a;) : 

*(a;) = exp {-iex^} {x±) , (^±1, (O, x\ x^) . (4) 

The stationary Dirac equation in the both representations is: 

(trPx + Ma^) 4'\x±) = etj4'\x^), = (0,Pi,P2) , (5) 
(a VP^a^ + Ma^) V^^"^) (x^) = eip^'^^ (x^) . (6) 

We note that the energy eigenvalues can be positive, £ = +£ > , or negative, e — < 0. 
We also can see that (5) and (6) are related by 

4-'\x±) = ^'^^-kx±) ■ (7) 

In what follows, we use the representation defined by C = 1- 

As the total angular momentum operator, we choose J = —id^ + (7^/2 that is the 
dimensional reduction of the operator in 3 + 1 dimensions. The operator J commutes 
with the Hamiltonian H. Therefore, we can consider Eq. (5) separately in each of 
eigenspaces of the operator J, 

Ji:^ {x^) ^(l-l^-^^ Vi'^ {xa) (8) 



Representing the spinors ipi^^ in the form 



^i^^ {x±) = 9i{^)ipi (r) , gi{^) = exp 



(9) 
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we reduce equation 5 to the radial Dirac equation for the radial spinor -0; (r) , 



n 



a" 



i<dr-\ 



l^ + l-l{l-a')+A{r) 



a 



(10) 
(11) 



with the radial Hamiltonian^ h; U defines the action of the spin projection operator on 
the radial spinor in the subspace with a given I, 

(rP±gi{(p)^i (r) = gi{(p)Ilil:i (r) . 

It is convenient to represent the radial spinor as 

^;(r)= [a'{e-U) + M]ui{r), (12) 

where 



uiir 



Vl 



(T=±l 



(13) 



and Ca are some constants. It follows from (10) that H^u — {e^ — M^) u; therefore, the 
radial functions 0/,ct('") satisfy the equation 



(f d p I 



UJ 
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e(/^+^-2(i-^; 



Ap 



>W = 0, (14) 



p = 7rV2, 7 = e|5|, { = sgnS, i/ = /x + / - (1 + a) /2, a; = £^ 



Solutions of the equation (14) were studied in [4]. The results can be summarized as 
follows. 

For any /, there exists a set of solutions (f)i^a = i<Pm;i,a- , m = 0, 1, 2, . . .), 



(t>m,l,a{f) = Im+\v\,m (p) , 



(15) 



square integrable and regular^ at r = 0. Here In,mip) ^-re the Laguerre functions that are 
presented in Appendix A. 

For / = 0, there exist solutions square integrable and irregular at r — ii p, ^ 0. A 
general irregular solution for / = and p ^ is 



a^p-il + a)/2, 2A = u;/7-e[/x-(l-a)/2] , 



(16) 



where are the Whittaker functions (see [38] , 9.220.4). The spinors in (10) that 

are constructed using these functions are square integrable for arbitrary complex A. The 

^By the radial Hamiltonian h, we mean the whole family of h with different I = 0,±1, ... and arbitrary 

B. 

^Here, we use the terms "regular" and "irregular" at r = in the following sense. We call a function 
to be regular if it behaves as r'^ at r = with c > 0, and irregular if c < 0. We call a spinor to be regular 
when all its components are regular, and irregular when at least one of its components is irregular. 
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functions t/jx ^ were studied in detail in [4], some important relations for these functions 
are presented in Appendix A. We see that interpretation of cu as energy is impossible for 
complex A. For real A there exist a set of solutions (16) which can be expressed in terms 
of the Laguerre functions with integer indices: 

Cl(^) = Im+,,-l,m (p) , (7 = 1, m = 0, 1, 2, . . . , 

C-iW =^m-^,™(p) , a = -l, m = 0,1,2,... . (17) 

All the corresponding solutions ipi{r) of Eq. (10) are square intcgrable on the half-line 
with the measure rdr. The Laguerre functions in Eqs. (15), (17) are expressed via the 
Laguerre polynomials. 

Eigenvalues u and the form of the spinors depend on sgnB. In what follows, we present 
the results for B > 0. The results for B < cannot be obtained trivially from the ones 
for B > 0. We present them in Appendix B. The spectrum of uu corresponding to the 
functions (pni,i,a{r') is 

/ 27(m + / + /x), /-(l + a)/2>0 , . 

\ 27(m + (l+a)/2), Z-(l + a)/2<0 ' ^^^> 

and the spectrum of u) corresponding to the functions ^.(r) is 

r 27(m + /.), a=l _ 

[ 27m, a — —1 ^ ' 

We require that the spinors ui (r) be eigenvectors for H, such that the functions Um^i^± 
satisfy the equation 

nMm,/,±(r) = ±^/ujum,l,±{r) . (20) 

We now can specify the coefficients in (13). 
In the case a; = 0, we have 

"»■'«=( L.M )•' ^ "»('■'=( °t A) )'' = °- 

This can be easily seen from the relations (118) - (121) for the Laguerre functions In,m{p) 
In the case a; 7^ 0, we have 

0m,O,l(^) 



m+l, 



,(r) ) . ' = 0. '^ = 27(m+ 1) ■ 



<±M=U*U-.W j'' = °.- = 2T("' + '')- (22) 

For 7^ 0, we construct solutions of the Dirac equation using the spinors u correspond- 
ing to the positive eigenvalues of the operator 11. These solutions are 

i^m,iir) = N [a^ {e-V^) + M] M„,,i,+(r), 1^0, 

i^'rA^ir) = N [a' {e-V^)+ M] z.^5(r), 1^0, (23) 



where is a normalization constant. Substituting (23) into (10), we obtain the two 
types of states corresponding to particles, ^ip, and antiparticles, with s = ±s ~ 
±\/M^~+1v, respectively. The particle and antiparticle spectra are symmetric, that is, 
|+£| = for the given quantum numbers m and I. 

We now consider the case u ~ 0. As it follows from (10) and (21), only the negative 
energy solutions (antiparticles) are possible. They coincide with the corresponding spinors 
M up to a normalization constant, 

-^o,i{r) = Nuo,i{r): I < -1; -V'o W = Mir): 1^0- (24) 

Thus, only antiparticles have the rest energy level. The particle lower energy level for 
Z < is +e = + 27. 

All the radial spinors '4'm,i{f) are orthogonal for different m. The same is true for both 
the spinors ip^^. In the general case, the spinors of the different types are not 

orthogonal. Using Eq. (124) of Appendix A, we can prove this fact and also calculate the 
normalization factor which has the same form for all types of the spinors. 



N = 




^ (25) 



e - ./uf + M2 



In addition, on the subspace / = 0, there are solutions of Eq. (10) that are expressed 
via the functions ipx,a{p) (16)- Wc represent these solutions as 

Mr)= [a'i6-U) + M]u^{r), 

Uoj{r) = ciu^^i{r) + c-iu^-i{r), Uu,,a{r) = ^uj,a{r)v,j ■ (26) 
Using the relations (130) for the functions ■?/'A,a(p), we obtain the useful expressions 

, LO 

nMu;,i(r) = iy/2'-fu^_x{r), Iiu^_i{r) = -i—=u^^i{r) , (27) 

V27 



Using Eq. (132) from Appendix A, we can see that the spinors ■0a;(r) and ■0w'(r), cu ^ cu', 
are not orthogonal in the general case. 

The completeness of all the obtained solutions is also an open problem. 

We conventionally treat the obtained solutions and spectrum as the respective energy 
eigenvectors and energy eigenvalues for the Dirac Hamiltonian H. In view of the confining 
uniform magnetic field and the apparent self-adjointness of H, we expect bound states 
and real discrete spectrum. What concerns the subspaces / 7^ 0, our expectations are 
completely realized. But in the subspace / = for fi ^ 0, there exist square-integrable 
solutions with complex eigenvalues (in what follows, we call the subspace / = and the 
subspace I ^ the respective critical and noncritical subspaces) . This implies that there 
is the problem of self-adjointness for the Dirac Hamiltonian with the magnetic-solenoid 
field, at least in the critical subspace. We solve this problem in the subsequent Sec. 3. 



2.2 Solutions in 3+1 dimensions 

To take the symmetry of the problem under the ^-translations into account, we use the 
following representations for 7-matrices (see [18]), 



„o 



\ 1 _ icr \ 2 _ / -icr \ 3 _ / 7 



^ = -<73 ' ^ = -ia^ ' ^ = ' ^ V -/ 
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In 3 + 1 dimensions, a complete set of commuting operators can be chosen as follows 

(7^ = -i7°7^7^7^): 



H, = -id3, = -id^ + EV2, = 7^7^ (M + -f^P^) /M . 
We require that the wave function be an eigenvector for these operators, 

= e^, 

53^ = sM/M<lf. 



(28) 



(29) 
(30) 
(31) 
(32) 



Here, M = s/WTl^, is the z-component of the momentum, and is the z- 
component of the total angular momentum. We note that the energy eigenvalues can be 
positive, e £ > , or negative, e = _e < 0. The eigenvalues j"^ are half-integer, it is 
convenient to use the representation = (/ — /q — |), where / = 0, ±1, ±2, .... To specify 
the spin degree of freedom, we choose the operator that is the 2;-component of the 
polarization pseudovector [42] , 



^ "2^ + ' ^ 2^ + ' 



(33) 



the eigenvalues of the corresponding spin projections are sM/M, s = ±1. 

Then, in 3 + 1 dimensions, we can separate the spin and coordinate variables and obtain 
the following representation for the spinors ^: 



*(x) = exp {-isx^ + ip^x^] "^six^) 



1 + Ip^ + sM) /M 



1+ (p' + sM)/M ^sAx±) 



(34) 



Here, ipe,s{x±) are two-component spinors, x± — {0,x^,x'^,0), and is a normalization 
factor. 

The equation (29) is then reduced to the equation 



(<tPx + sMa^^ ^e,s{^±) = eij,,s{xi_), Pi. = (0, Pi, P2, 0) 
Representing -^^^^ [x±) in the form 

where gi {ip) is given by Eq. (9), we obtain the radial equation 



(35) 



(36) 



(37) 



where hs is the radial Hamiltonian acting on the subspace with the spin quantum number 
s, n is given by Eq. (11). We note that 



'4^e-l{x±) = Cr^1p-e,l{x±) 



(38) 
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Wc can sec that for fixed s and p'^, Eq. (35) is similar to Eq. (5) in 2 + 1 dimensions. 
Therefore, after separating the angular variable by (9), the radial spinor (r) (36) can 
be obtained from the radial spinor ipi (r) (9) by substituting M by M. The same is true 
for the particular case I — 0. Here, the radial spinor {t^) can be obtained from the 
radial spinor ■0^ (r) (26). 

Using the results for the (2 + l)-dimensional case, we conclude that in the critical 
subspace, the complex eigenvalues of Eq. (29) do exist if /i 7^ 0. This means that the 
abovementioned problem of self-adjointness for the Dirac Hamiltonian in 2 + 1 dimensions 
is reproduced in 3 + 1 dimensions. 

3 Self-adjoint extensions 

We here solve this problem using von Neumann's theory of the self-adjoint extensions of 
symmetric operators [43]. The canonical procedure includes the following steps: initially 
defining the Dirac Hamiltonian i7 as a symmetric operator, then evaluating its adjoint 
and its closure H, then finding the deficiency subspaces D"*" and T>~ and the deficiency 
indices, and, at last, in the case of equal deficiency indices, describing the isometries from 
T>'^ to V~ that define the self-adjoint extensions of H. This constitutes von Neumann's 
theory that is applicable to the general case. In Sec. 5, we show that in our particular 
case, this procedure can be significantly reduced: it is sufficient to evaluate W and then 
analyze its symmetry properties; of course, we thus follow one of the main ideas of von 
Neumann. 

3.1 Extensions in 2+1 dimensions 

We first study the (2 + l)-dimensional case. This case was partially considered in [37]. 
The results of [37] are recovered by our consideration based on the detailed study of the 
Dirac equation solutions. We also generalize these results for an arbitrary sign of B, 
which allows determining the nontrivial dependence of the spectrum on the signs of B 
and $. The obtained results are then extended to the (3 -|- l)-dimensional case. From 
the mathematical rigor standpoint, our exposition is rather qualitative, some important 
technical details can be found in Sec. 5. 

The problem is to define the Hamiltonian (1) as a self-adjoint operator in the Hilbert 
space of square-integrable two-spinors ip{x±). By separating variables, Eq. (9), the 
problem is reduced to the corresponding problem for the radial Hamiltonian h (10), (11) 
for every I — 0, ±1, ... in the Hilbert space of two-spinors ip (r) square integrable on the 
half-fine with the measure rdr. We start with the choice of the initial domain T> (h) for the 
operator h. Because the source of the problem under consideration is the singularity of 
the AB potential at the origin r = 0, we first try to avoid the troubles associated with this 
singularity. Therefore, let V {]%) be the (sub)space of absolutely continuous two spinors 
'0(r) vanishing sufficiently fast as r ^ 0; of course, hiplr) must be square integrable 
together with ip{r). It is easy to verify that the initial h is symmetric by integrating by 
parts. Its adjoint is given by the same expression (10), (11), but is defined on the 
domain T>{h^) D T>{h) of absolutely continuous two-spinors not necessarily vanishing as 
r ^ in the critical subspace / = for 7^ 0. Its closure h is defined on the domain 
T>{h) of absolutely continuous two-spinors vanishing as r — > 0. 
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Then we have to find the deficiency subspaces and , =Ker(/i"l' =]= iM) (here, 
M is introduced by dimensional reasons), and the deficiency indices n± {h) = dim(P=''), 
i.e. to find the number of hnearly independent square-integrable solutions of the equations 

/it^±(^) ^ ±iMi/j^{r), = W + a^M , (39) 
n^^-i!^dr + y fx + l-^{l-a^)+A{r) (40) 

In the noncritical subspaces / 7^ 0, there are no such solutions. In the critical subspace 
/ = 0, if /X 7^ 0, there is only one solution for (r) and ip~ (r), namely, 

^Hr) - n(^^%,.^,),B>0, (41) 



(42) 



where 

Mr) = V'A,a(p), 2A = -2MV7 - e (/X - (1 - (7) /2) , a = ±1 , 

see (16), and N is the normalization factor such that the functions ?/'^(r) are normalized 
to unity. If = 0, there are no such solutions. This means that in the non-critical 
subspaces, the deficiency indices are (0,0), whereas in the critical subspace, these are 
(1, 1) unless /X = 0; if /X = 0, the deficiency indices are (0, 0). Therefore, in the noncritical 
subspaces and, ii /i — 0, in the critical subspace as well, the radial Hamiltonian h is 
essentially self-adjoint, i.e., its unique self-adjoint extension is its closure h = . This 
guaranties the completeness of the corresponding solutions obtained in Sec. 2. In the 
critical subspace, the radial Hamiltonian with 7^ has a one-parameter family {h^} of 
self-adjoint extensions, which is homeomorphic to the group U{1). Each member of 
this family is defined by the isometry ■0+(r) — > e^^'^~{r) from onto V~ , < Q < 27r. 
The domain of hP is 

V{h^) = {x{r)=ij{r)+c[i)+ {r)+e'^ij- ir)]-. (r) e V(h)} , c e C , (43) 

We now note that V (/i^) is completely defined by the asymptotic boundary conditions 
on the two-spinors x (^) = ix^i^), x'^i^)) as r — >• 0. Namely, because -j/^ (r) — >• as r — 0, 
then, if c 7^ 0, the behavior of x (r) as r — > is defined by the singular behavior of i{j~^{r) 
and ip~{r) near the origin. Using the behavior (131) of the function ipx^a{p) a-t small p, 
we find 

xHr){Mr)' ^ ^ I (tan f -i)r(,)r(i+MV7) U j ' ^ > " 

[ (tanf-l)r(M)r(l-M+MV7) V T / ' 

if c 7^ 0. (Of course, c can be equal to zero, then x(r) — > as r — > 0.) These alternative 
possibilities are just the asymptotic boundary conditions that uniquely specify the domain 
V (/i^) and thus the self-adjoint extension h^. The self-adjoint asymptotic boundary 
conditions can be consolidated into a single formula 

X (r) = c '"^(i^p"' ) + 6 (r) , r ^ , (45) 
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where iA is given by the right hand side of (44), —oo < A < oo, e (r) ^ as r ^ 0, and 
c is an arbitrary constant. We can verify that in the hmit 7—^0, the right-hand sides of 
(44) coincide with the corresponding expressions obtained in [11] for the case of pure AB 
field. In fact, the family of self-adjoint asymptotic boundary conditions is the same as 
in pure AB-field case, and is independent from B. These conditions are dictated by the 
singular behavior of the AB potential at the origin. 

For our purposes it is convenient to pass from the parametrization by Q to the parametriza- 
tion by the angle 6, < © < 27r, such that 

lim = ^ tan -- + - (46 

if c ^ 0. 

Therefore, the solutions (26) obtained in Sec. 2 must be subjected to the asymptotic 
condition (46) as r 0, which guaranties the orthogonality and completeness of the 
corresponding solutions. Using (27), (139), and (131), we find 

tan ( ^ + = ^ 2.r(.)r(i-./2.) [ , ) ^ ^ > " _ 

V 4 2 / 1 M r(i-^)r(i-c^/27) f mI\ ' ^ < 0. 



(s-M) 2A'-lr(/i)r(l-/i-a)/27) 7 

3.2 Extensions in 3+1 dimensions 

We now pass to the (3 + l)-dimcnsional case. The hclicity operator Sh = SP/|P| is 
commonly used as the spin operator. It is related to the zero-component of the polarization 
pseudovector (33) by Sh — S^M/\'P\. Extensions, in which the operator commutes 
with the Hamiltonian, were constructed for the particular case = [24, 25] . However, 
the set of such extensions does not exhaust all the possible Hamiltonian extensions in 2 + 1 
dimensions. Our choice of the operator as the spin operator allows us to separate the 
spin variables from the beginning and to remain with the extension problem for the radial 
Hamiltonian only. Therefore, we can efficiently apply our experience in 2 -|- 1 dimensions 
to the (3 + l)-dimensional case. 

In the (3 + l)-dimensional case, after separating the ^-variable, the Hamiltonian H 
must be defined as the self-adjoint operator in the Hilbcrt space Ti of four-spinors of the 
form (34). This Hilbert space can be represented as the direct sum of two orthogonal sub- 
spaces labelled by the spin quantum number s: H = © {^'-i}- These subspaces are 
invariant with respect to H, and consequently, the Hamiltonian (1) can be independently 
considered in each of the subspaces. Using Eqs. (35) and (9) allows reducing the prob- 
lem to the radial Hamiltonians hs (37) acting on the subspaces of the two-spinors with 
the given spin quantum number s = ±1,... and orbital quantum number / = 0, ±1,.... 
The problem with the radial Hamiltonian /i^in the (3 -|- 1) -dimensional case is absolutely 
similar to the problem with the radial Hamiltonian h in the (2 -|- 1) -dimensional case and 
is solved by the same von Neumann's method. Because the procedure for defining hg as 
the self-adjoint operator literally repeats the procedure for h in the previous section, we 
only outline the main steps. 

We first define the radial Hamiltonian hs as a symmetric operator choosing for the 
initial domain T> (hg) the space of absolutely continuous two-spinors vanishing at the 
origin. Its adjoint hi and its closure hg are described just in the same terms as the 
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respective W and h from the previous subsection. We now apply von Neumann's theory 
to each of the subspaces. To find the deficiency subspaces and deficiency indices of the 
operators hs, we have to solve the equations 

h\il;f (r) = iisMV^f (r) , h\^Ii^ + sMa\ s = ±1 , (48) 

where 11^ is given by (40). These equations are the copies of Eq. (39). Using Eqs. (41), 
(42) and (38), we find that for Z = and 7^ 0, the solutions are 



'Pt{r) = N[ ).B>0, (49) 



M 



0s,+i(r) 



«fM = iV| ±,,.w.M^,_,M ].B<0, (50) 

0s,a(?^) = ^X,a{p) ' a = M - (1 + 0-) /2 , 

2A = -2MV7-e(/x-(l-<j)/2) , a = ±l, 

whereas for Z 7^ and for Z = and = 0, there are no square integrable solutions, which 
means that for each s = ±1, the deficiency indices in the noncritical subspaces I ^ are 
(0, 0), whereas in the critical subspace / = 0, the deficiency indices are (1, 1) unless = 0; 
if = these arc (0,0). This implies that in the noncritical subspaces and, if = 0, 
in the critical subspace as well, the radial Hamiltonian hs is essentially self-adjoint, i.e., 
its unique self-adjoint extension is its closure h — h\ whereas in the critical subspace, 
if // 7^ 0, there exists a one-parameter family {h^"} of the self-adjoint extensions of the 
radial Hamiltonian hg, labelled by the parameter , < fig < 27r. The domain of h^" is 

^ {h^') = {Xs (r) = i's (r) + c [^+ (r) + e^^^C (0] : (^) e V (hs) } , c G C . (51) 

Quite similarly to the (2 + l)-dimensional case, h^" is specified by the self- adjoint 
asymptotic boundary conditions on the two-spinors x{f^) = (X"'^(^)i X^(^)) as r ^ 0. With 
the parametrization by the angle 0^ similar to (46), these boundary conditions are 

xMr)(Mr)''' 0X 

lim = s« tan - + — ^ , s = ±1 (52) 

xUr) [Mr) V4 2 / 

if c 7^ and Xs(r)— >Oasr— >Oifc = 0, which can be consolidated into a single formula 

iAs (MrY 

Xs{r) = Cs\ ^\ ^L^ |+es(r),r^O, 

Mr 



where iAg is given by the right hand side of (52), —00 < Ag < 00, e (r) — > as r ^ 0, 
and Cs are arbitrary constants. 

Therefore, in each subspace s = ±1 the above solutions '^l]^^^s ir) in the critical subspace 
I — must be subjected to the condition (52). 

The remarks on the orthogonality and completeness of the corresponding solutions are 
similar to those in previous Sec. 3.1. 
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The final result is that in 3 + 1 dimensions, there is the two-parameter family of the self- 
adjoint Dirac Hamiltonians. This family is the manifold U{1) xU{l). We emphasize that 
this is true only if we require that be conserved. If we don't require the conservation of 
S^, then it follows from the above consideration (for example, any vector satisfying Eqs. 
(29)-(31) can be constructed using two orthogonal vectors determined as a superposition 
of ^^+1 and \E'_i) that the deficiency indices in the critical subspace are (2, 2) unless ^ = 0. 
Consequently, the manifold of self-adjoint Dirac Hamiltonians is the group U (2) and is 
four-parametric. Only for the submanifold U{1) x U{1) C U{2), the Hamiltonian H and 
the operator have a common set of eigenf unctions. 



3.3 Spectra of self-adjoint extensions 

We now study the spectra of the self-adjoint extensions h^. To find these, we have to solve 
the transcendental equations (47) for cu considering two branches of e, one for particles 
and another one for antiparticles, ±e — ±y/ + u. Introducing the notation 

cj = 27a;, x= ,x={,e^- M^)/2^,Q{x) = ^ + = ±, 
for S > 0, we can rewrite Eq. (47) as 

Given a; for 5 > 0, we can obtain a; for 5 < with the substitutions 

fj (^) ^l/fj{n), ix^l- n. 

In what follows, we therefore consider the case B > Q only. 

The possible solutions x = x{r]) of the equation (54) are the functions of the parameter 
77 (i.e., of ji, 7/M^, and ©) and are labelled by m = 0, 1, ... . We find the following 
asymptotic representations for these solutions as |?7| — > 0: 

/ \ A A ('''"A'') + 1 — /i) 

Xm [r])^m + iXXm, i^Xm = ^, . ?7 , m = 1, 2, 3, . . . , 

7rr(m)(y (m) 

All Xm (0) , m = 1, 2, ... are positive and integer. The asymptotic representation of ^xq {rj) 
as 1 77 1 — > is discussed below. The function ^xq {vj) vanishes at the point rj = 2T{n) , and 
in the neighborhood of this point, it has the form 

/ N r(/x) - rj/2 

Here, iIj{x) is the logarithmic derivative of the gamma function T{x), and —'(/'(I) — 0.577 
is the Euler-Mascheroni constant [44]. As I77I —>■ oo,we find the following asymptotic 
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representations: 



rXr, 



(^77) ^m + ii + Axm, m = 0, 1, 2, . . . , 77 

(^■77) = ??7, - 1 + /i + AXrn, 772 = 1, 2, 3, . . . 

sin (TT/y) T{m + (tti + /i) 



00 



?7 



-00 



(57) 



These approximations hold only for |Aa;m| <^ yuand \xo ^ /i. 

According to [45] (see Corollary 1 of Theorem 8.19 therein), if Ti and T2 are two self- 
adjoint extensions of the same symmetric operator with equal finite deficiency indices 
{d,d), then any interval {a,b) C TZ not intersecting the spectrum of Ti contains only 
isolated eigenvalues of the operator T2 with total multiplicity at most d. Wc take the 
extension with O = 7r/2 whose eigenvalues are = Ma/ 1 + 27 -^xq (00) /M^ and 
= ±M-\/l + 2^ ±Xm (±cxo) /M^, m > 1. The above theorem implies that if (a, 6) is 
an open interval, where a and 6 are two subsequent eigenvalues of h^^ with = tt/2, 
or ±£ = 0, then any self-adjoint extension /i^ at © 7^ n/2 has at most one eigenvalue 
in (a, 6). According to [46] (see Theorem 3 from Sec. 105 of Chapter VIII therein ), 
for any e G {a,b), there exists a self-adjoint extension with the eigenvalue e. As it 
follows from (54) and (57), in the ranges {m — 1 + n < ±Xm {rj) < m + fi, m > 1) and 
(-MV27 <+ xo{ri) < n), the functions ±x{ri) = {±e^ - M^) /2-f are one-valued and 
continuous. This observation is in complete agreement with the above general Theorems. 
The functions ±Xm (t)) were found numerically in the weak field, 'j/M'^ <^ 1, for some 
first TTi's. The plots of these functions (for // = 0.8) are presented in Figs. 1 and 2. 



x=co/2y 




Figure 1: Particle lower energy levels in dependence on the parameter 77+ = 

r(l-M) \2M'^) ^'^^^ 14^2^ 

We can see that 5xm — x^+i (?7) — {v) 1 with increasing m. It follows from the 
equation (54) that 



SXm - 1 = TT ^ {cot {iTXm) - COt [tT {x^ "A*)]} ^ 



1-// 



- (5g , 771 > 1 , (58) 
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x=co/2y 





Figure 2: Antiparticle lower energy levels in dependence on the parameter rj- 
r(i-M) \2M2) ^"^^ U ^ 2; 



where SQ = ^lnQ(x)|^^^ < l/x^- The curve X5 (77) can give an idea about the 
behavior of the functions (v) for large m . 
In the sequel, we discuss some limiting cases. 

We first consider the weak fields B for which ^ 1 and the nonrelativistic elec- 

tron energies, a;^ (?]) ^ 1. Here, the functions ±x (rj) change significantly in the 

neighborhood of 77 = only. Beyond the neighborhood of 77 = 0, the functions ±x {rj) 
take the values close to the corresponding asymptotic values given by (57). 

In the ultrarelativistic case, Xm{Tl)'~f /M^ ^ 1, the behavior of Xmifj) qualitatively 
depends on /i . One can distinguish the three cases: < /i < 1/2, /i > 1/2, and n = 1/2. 
If < /i < 1/2, then the interval near 77 = on which the functions change significantly 
narrows down with increasing m. If // > 1/2, then this interval expands with increasing 
777. . For 11 — 1/2 and — 1<(| + ^)<|, we find the asymptotic representation 



<;Xm (^) = "I + ^ ( ^ + ^ ) , 771 > 1 . 



27r/ 



(59) 



We see that negative ±Xq {rf) exist only for 7/ > 0. This means that in the problem 
under consideration, there exist only one particle state and only one antiparticle state 
with energies \e\ < Mfor7r/2 < < 37r/2. The same situation was observed in the pure 
AB field case [11]. The minimal admissible negative Xo{r)) is defined by the condition 
e = 0. For the strong fields B for which 7/M^ ~ 1, the quantity Xq (7/) is close to zero. 
Let Oo correspond to such an extension that admits £ = 0. The value of 60 is defined by 



tan 



TT Go 

4 2 



r(i-/x)r(/x + MV27) 
2'^r(Ai)r(i + M727) 



M2 

7 



i-n 



(60) 



In the weak fields, 7/M^ <^ 1, xq (77) takes large absolute values, and the angle 0o is 
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defined by 

and independent from tlie magnetic field. It follows from (60) that in the superstrong 
fields B for which ^ 1, the angle ©0 is also independent from the magnetic field. 

In the weak magnetic fields, <^ 1, and for the nonrelativistic energies, xq^/M"^ <^ 

1, we obtain the relations 

+a;o (r/) = - , (62) 

-Xo{r,) = -{r,M^hy'' (63) 

that are valid if rj in (62) is small and r^M"^ 3> 1 in (63). 

We now consider the particular case G = — 7r/2. It follows from (47) that for > 0, 
there exists _e = —M. The energies \e\ > M are defined by the poles of r(l — x) or 
r(l — fj, — x) ior B > or B < 0, respectively. The spectrum e coincides with the one 
defined by Eqs. (137) and (23) for ijj^. Moreover, using the relation (128), we can see 
that the spinors ■0w(r) coincide with up to a normalization constant, 

il)^ (r) oc (r) for = -tt/2 . (64) 

In the case O = 7r/2, we have the following picture: It follows from (47) that for B < 0, 
there exists = M. The energies |£| > M are defined by the poles of r(/i — a;) or r(l — x) 
for > or i? < 0, respectively. The spectrum e coincides with the one given by Eqs. 
(137) and (23) for t/j^^. It follows from (128) that the spinor ■0w(r) coincides with ip^^ up 
to a normalization constant, 

^^(r) oc (r) for 9 = 7r/2 . (65) 

Using the results for S < 0, which are presented in Appendix B, we can conclude that 
the spectrum asymmetry holds for the spinning particles in the magnetic-solenoid field. 
There is a relation between the three-dimensional chiral anomaly and the fermion zero 
modes in a uniform magnetic field [39] (for review, see [40, 41]). We see that the effect 
also holds in the presence of the AB potential. 

The spectrum asymmetry is known in 2-1-1 QED with the uniform magnetic field. In 
the uniform magnetic field, the states with u = for / 7^ are observed if sgn/ = —sgnB 
(for the antiparticle ii B > and for the particle ii B < 0). The spectrum changes 
mirror-like with the change of the sign of the magnetic field. We see that for I ^ 0, the 
spectrum properties in the magnetic-solenoid field is similar to the spectrum properties 
in the uniform magnetic field. The presence of the AB potential is especially essential for 
the states with / = where the particle penetrates the solenoid. 

The spectra in 3 + 1 dimensions can be obtained from the results in 2 + 1 dimensions. 
Namely, we use the fact that the solutions 'ipe,i{x±) in 3 + 1 dimensions arc obtained from 
the solutions ijji^\x±) in 2 -|- 1 dimensions. Therefore, the spectra in 3 -|- 1 dimensions 
are obtained from the results in 2 -|- 1 dimensions by substituting M by M and using the 
relation (38). As a consequence, we obtain an additional interpretation of Figs. 1 and 2. 
In particular, Fig.l shows the lower energy levels for particles with spin s — 1, and Fig. 
2 shows the lower energy levels for particles with spin s = — 1. 
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4 Solenoid regular izat ion 



4.1 Spinning particle case 

We can introduce the AB field as a limiting case of a finite-radius solenoid field (the 
regularized AB field), which allows fixing the extension parameters. This approach to the 
pure AB field was first proposed by Hagen [12]. In what follows, we consider the problem 
in the presence of a uniform magnetic field. For this, we have to study the solutions of the 
Dirac equation (1) with the combination of a finite-radius solenoid field and a coUinear 
uniform magnetic field. 

Let the solenoid have a radius R. We assume that inside the solenoid, there is an 
axially symmetric magnetic field B'^'^{r) that creates the flux ^ — [Iq + jj) $o, $o = 27r/e, 
such that e /^f -B*" (r)rdr — Iq + ii, and outside the solenoid (r > R), the field B^^(r) 
vanishes.. The function B'^^{r) is arbitrary but such that the integrals in the functions 
d {x) and h [x) in (71) are convergent. We take the potentials of the field B^"'{r) in the 
form 

eAT = ^ (x) eA, = (x) ^ , (66) 

where 

rx 

■d{x)^ I f {x') x'dx', f{x) = R'^eB'''{xR), x^r/R. 
Jo 

The potentials of the uniform magnetic field are 

Ao^O,A,^A (r) A, = -A (r) ^ , A (r) = Sr72 . (67) 

Outside the solenoid, the potentials have form (3). 

We analyze the solutions of the Dirac equation in the field defined above . For this, 
we have to solve the equation inside and outside the solenoid and continuously sew the 
corresponding solutions. We call the corresponding Dirac spinors the respective inside 
and outside solutions. 

We first study the problem in 2 + 1 dimensions. We require that the solutions be square 
integrable and regular as r — 0. In the same way as in Sec. 2, we obtain that the inside 
radial spinors il^'^^^iir) (r < R) satisfy the equation 



h'^'i^Z (r) = (r) , h}'' = n- + a^M , 



where 



3 



R X 



l-lo-\{l-'J^)+'&{x)+^PRX^ 



I a\ PR = 7i?V2 . (68) 



We require that the functions ip^^i{r) be square integrable on the interval {0,R)- For 
LU = (1^1 = M), the solutions are 

+V'S(r) = 05,1(^)^^1, l-lo>l, 
-V'S(r) = 05,_i(x)^;_i , /-/o<0, 

<Po,iA^) = ca;!''! exp i^a x'^ {'d{x) + Cprx^) dx^ t) = I - k - {1 + a) /2 , (69) 
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where c is an arbitrary constant. For a; 7^ 0, we represent the spinors in the form 

where c^- are arbitrary constants. The functions 0^'^ {x) satisfy the equation 

<IC {x) = (70) 



1 ^ O -1 



and must be regular at r = in order to satisfy the square integrabihty condition for 
'^oj'i ('")■ prime interest is in the hmiting case R ^ 0. For our purposes, it is enough 
to use the approximation pji <^ 1, ujR^ ^ 1. Rejecting the terms proportional to in 
(68) and (70), we find that solutions of Eq. (70) are 

,in( N _ / cx\'^\e"^^''\ ar]>0, 

b{x) = f dxx~^^{x) . (71) 
Jo 

The outside solutions (r > R) satisfy the equation 

^C? (r) = er.1 (r) (72) 

and must be square integrable on the interval {R, oo). Here, h is defined by Eqs. (10) 
and (11). The general form of the outside solutions is 

(^) = K - n) + M] {c^cpif{r)v^ + tc_^cP1%{r)v_,) , 
C*W=^A,a(p) ,a = l + ^l-{l + a)|2,2\ = uJh-i{l + ^i-{l-a)/2) (73) 

The solutions -0°"/ (r) and V'^"; (r) must be sewed continuously at r = i?, 

[R) = [R) , (74) 

and must satisfy the normalization condition 

j^in , j^out 
pR POO 

= / (V'wW)Vw(^H^, / {^:1ir)y^::l{r)rdr. (75) 

JO JR 

We can treat the AB field as a limiting case of the finite-radius solenoid field if 

limTV^^O. (76) 

We can realize the sewing condition (74) imposing the following conditions on the functions 

(''") ^i^a i"^) at r = i?: 

(i? - e) = (i? + e) , ^0 (i? - e) = -^0 (i? + e) . (77) 

dr dr 
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It is convenient to use the representation (128) for -0;^ ^ (p) in (73). The functions (pf'^ (r) 
then are 

1 — CX 1 ~t~ 
<l>r^ il') = 0^4. ,m. (P) + KIma,n. (p) , = A m,, = A — , 

a^r — K smno-TT, — —K smmo-TT, K — — , (78) 

sin (rio- — 777,0-) 

where n„, nia- are real numbers. 

Using (77), we can find the coefficients a^., b,: for the case I — Iq <0, 

ai = Pn^'^'-'^^'cd,, h = pI^'^'^'^cI, , (79) 
a-i = p-^^'-^/'^'ca.,, 6_i = p(;+-)/^cLi , (80) 

whereas for the case I — Iq > 0, 

ai = Pn^'^'-'^^'c%, b, = pr'^-^^/^+^c'ft; , (81) 
a-i = p;i^'-^''^/'c'a'_„ = pg+-)/^c'6'_, , (82) 

where the non-vanishing coefficients a, b, a', b' are independent from and the coefficients 
c and c' are normahzation factors which depend on p^. 

Calculating the normalization factors, in the limit i? ^ 0, we obtain 

Oi = const 7^ 0, bi = 0, a_i = 6_i = 0, / > 1 , 

oi = 0, bi — const 7^ 0, a_i = 6_i = 0, / < 

for / — /o < and 

oi = 6i = 0, a_i = const 7^ 0, 6_i = 0, / > , 

ai = bi — 0, a_i = 0, 6_i = const 7^ 0, I < —1 

for Z — Zo > 0. For / = 0, the value of the coefficients is defined by sgn$. We can verify 
that the condition (76) is satisfied. 

We thus obtain that for any sign of B, the solutions are expressed via Laguerre poly- 
nomials (23). In particular, for / = 0, we find that the solutions (r) coincide with 
either (r) or (r) in accordance with sgn ($), 

Lout /X _ f i^L (r) , sgn ($) = +1 



(r) - I ^7} ' ($) = _! ■ (S3) 

In Sec. 3, we have found the relation between the extension parameter values and the 
types of solutions in the critical subspace I — (64), (65). We are now in a position 
to refine this relation. Namely, if we introduce the AB field as the field of the zero- 
radius limit for the finite-radius solenoid , then the extension parameter is fixed to be 
= —sgn ($) tt/2. In addition, this way of introducing the AB field implies no additional 
interaction in the solenoid core. 

To solve the problem in 3 -|- 1 dimensions, we use the results in 2 -|- 1 dimensions 
presented above. In the limit — > 0, the solutions in the critical subspace are 



*r(3^±) = ^ V ~\ ^1 h 84) 

-1 + + sMj /mJ ^0 (<^) W 
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where the functions go{^p), ip^^f (r) arc defined in (9) and (83). (23). respectively. The 
values of the extension parameters in 3 + 1 dimensions are specified to be 




The interpretation of other possible O's via the limiting process for other regularized 
potentials is not reached so far. 

4.2 Spinless particle case 

For completeness, let us consider the regularization problem in the spinless particle case. 
For this case the Klein-Gordon equation with the magnetic-solenoid field is, in fact, re- 
duced to the eigenvalue problem for the nonrelativistic two-dimensional Hamiltonian. 
Therefore, the self-adjoint extension problem as well as the solenoid regularization prob- 
lem is similar for the relativistic and nonrelativistic case. 

From the classical paper [35] on, the AB effect in the spinless case was always associated 
with the radial functions regular at r = 0. However, the linkage between these kind of 
boundary conditions and a regularizing solenoid is an interesting important problem. The 
aim of this subsection is just to study this problem. 

One ought to say sclf-adjoint extensions of the nonrelativistic spinless Hamiltonian with 
the AB field were studied in several articles [6, 7, 8, 9]. The case of the magnetic-solenoid 
field was considered in [34] where the most general four-parameter family of admissible 
boundary conditions was obtained. The solenoid regularization with some particular 
distributions of the magnetic field inside the solenoid was studied in [6, 12]. As we know, 
the regularization problem with the arbitrary field inside the solenoid was not solved. 

Unlike the Dirac equation, the Klein-Gordon equation was not solved explicitly for the 
arbitrary field inside the regularizing solenoid. However, one can obtain some properties 
of the corresponding solutions without the explicit solving the Klein-Gordon equation. 
Thus, we can demonstrate that for the arbitrary field inside the solenoid, the hfting of 
the regularization {R 0) corresponds to regular radial functions only. 

To do that we use the above formulated in the Dirac case method. Solutions of the 
Klein-Gordon equation with a given energy e and angular momentum / — Iq have the form 



We find the radial function (r) sewing solutions inside and outside the solenoid. 
The outside solutions have the form [4], 



where, as before, u; — e"^ — M^. The critical subspace is defined by Z = 0, — 1. The inside 
solutions satisfy the equation (70) where one has to set a = 0. The approximation pn <^ 1 
and ujR'^ ^ 1 can be applied in this case as well. Therefore, we arrive to the following 
equation for the inside solutions. 



(f){x) = A^e^"^°e^('-'°)<^0(r) . 



0°"*(r)=V'A,a/2 (p)=p- 



^^^Wx,a/2 (p) , 2A = a;/7 - Ca, a = / + A* , 



(86) 



(87) 
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where < a; < 1. We are looking for solutions 0*" {x) that are regular at x = 0. Applying 
the sewing conditions (77), we obtain 

a<f>in (1) + 9,0^" (1) = 2aap''/ , a<j>in (1) - 9,0^" (1) = 2a6p-"/' (88) 

in the lower order in pR. 

For our purposes it is important to demonstrate that the following condition holds true 

|a| 0*^1) +5.0'" (1)^0. (89) 

To this end, let us find solutions of equation (87) that are regular at x = 0. The function 
i} (x) is analytic on the interval < a; < 1 as obeying the conditions listed in Sec. 4.1. 
Then (87) is the homogeneous ordinary differential equation with the regular singular 
point X = 0. It is known from the general theory [36] that there exist solutions of (87) 
that can be represented in the form 

oo 

0'" (x) = xl"! Yl ^kX^ ' (90) 

fc=0 

where Oq is an arbitrary constant and the other coefficients are defined via recurrent 
relations. The series in (90) is absolutely and uniformly convergent for < x < 1. Let 
us suppose now for simplicity that / {x) does not change its sign for < x < 1. Then we 
represent 0*" {x) in the form 

^ - / ^''''e-l'^^V^-) (^) , sgn (/) 77 < 



where the functions 9?*^^^ (x) obey the equation 

dx^ ' ^ ' ' " ' ^ xdx 



i/(±VW = 0, //W = //f^-|/(x)| , i/i^) = -^ + (l:T:2|77|-2|79(x)|)--^ (92) 



and the conditions 

^(-)(0) = 1, a:-2|''I^W(a;)|^^^ = l. (93) 

It is convenient to introduce into the consideration the retarded propagator G^gj (x, y) 
obeying the equation 

Et^G^^^{x,y)^?>{x-y) . 

It can be represented as G')^^ (^) u) = ~ v) G^^^ {x, y) where the functions G^^^ obey 
the conditions 

G^^^{x,x)^Q, ^GW(x,y) =1, 

x=y 

and can be found explicitly in the case under consideration, 

(x, y) = f^i^:^ ' exp {2 \h (x) | - 2 \h {y) \ } dx . (94) 

The differential equation (92) with the boundary conditions (93) is equivalent to the 
following integral equation, 

{x) = ^"t^ {x) + r (x, y) 1/ (y) | {y) dy , (95) 
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where 

(^^-^ (x) = 1, (x) = r x^\^\-^e^\''^^^\dx. 

Jo 

The solutions of the equation (95) can be found by iterations, 

oo 
k=0 

(x) = (x, y) 1/ (y)| (y) dy, Fo^^^ (^) = ^ (^) ■ (96) 

Every member of the series (96) is positive. The series converges uniformly. Then, it 
implies <^*^^^ (x) > and dxip^"^^ (x) > for x > 0. Therefore, 

0- (1) > 0, 9,0- (1) = -sgn (/) (1) + e-l^(i)la,^W (1) . (97) 

If sgn(/) a < 0, then (1) > and the condition (89) is satisfied. If sgn(/) a > 0, 

then it follows from (97) that 

|a|(/.^"(l)+a,0-(l) >0, 

and the condition (89) is satisfied as well. The obtained result can be extended to the 
function / (x) alternating in sign. In this case the interval [0, 1] can be divided into 
subintervals on which / (x) does not change its sign. The solutions 0*" (x) can be found 
successively on each such a subinterval beginning from the point x — 0. 
Applying the normalization condition for the sewed solutions, 

Jo J pa 

we get 

6 = 0ifQ;>0; a = OifQ;<0 (98) 

in the limit i? — > 0. Thus, the only regular at r = solutions remain in the limit it! — > 0. 
The conditions (98) also define the spectrum. Finally, it follows from (78) that in the 
limit — > the radial functions and the related spectrum has the form 

4'm,i (r) = Im+\i+t^lm {p) , ^ = 1 C^m + \l + fx\ + ^ {I + fx) + 1) . (99) 

5 Reduced self-adjoint extension method 

We here show that the general self-adjoint extension method can be significantly reduced 
for the radial Hamiltonian h given by Eqs. (10), (11). 

Wc recall that the problem is to define the formal matrix differential operator (10), (11) 
as a self-adjoint operator in the Hilbert space of two-spinors x (^) square-integrable on 
the half-line = {r > 0} with the measure rdr. It is convenient to pass to the standard 
measure dr on with the substitution x ('") = r~^/^^ (r) . The radial Hamiltonian then 
becomes (we do not change the notation for the transformed Hamiltonian h) 

h = U + a^M,U=( ° """A^]: (100) 

y a-i,-B J 
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where v — ^ -\- 1 — 1/2 and the operator Oj^^b: L2 L2 is 

5r + - + J = -ia_y_Bdr(yy,B (101) 

with 

feBr'^\ 

a^^B = r exp I I , a_^_Bcru,B = 1 • 

This h (100) must be defined as a self-adjoint operator in the Hilbert space 

n = L2 iR+) © L2 iR+) 

of two-spinors 

*(^)=U^(r) 

square intcgrablc on /?+ . 

We first note that o^M is a bounded self-adjoint operator in Ti. Therefore, the problem 
is equivalent to the problem of defining the apparently self-adjoint^ operator 11 (100) as 
a really self-adjoint operator in 7i : we must ensure the equality H = by the proper 
choice of the domain 

V{Ii) = V{a_,,_B) © V{a,^B) C H , 

for any z/ and B. The problem is thus, reduced to the problem of properly defining a^^s 
(101) as an operator in L2 (-R+) , i.e., defining V^g^^b) C L2 (-R+) , together with a^^^^B ■ 

We canonically start with defining 11 as a symmetric operator in 7i by taking the initial 
domain to be I^(n) — D ® D, such that the both a,y^B and Q-^-b are initially defined on 
D G L2 {R+) , where D is the linear space of C°°-functions with compact support. We 
note that 11 is densely defined because D <Z L2 {R+) and symmetric, 11 C H^, which is 
simply verified by integration by parts. 

The next step is the evaluation of its adjoint . It is evident that 

= (| j , l^i nt) = V{al^) © V{al^^_^) , (102) 

for any u and B. Because Gi^^b is a simple differential operator in L2 (-R-i-) , its adjoint a|,g 
is evaluated simply (by the standard method for differential operators in L2 {R+))- its 
form is 

4,6 = -icru,B dr(T-i,-B (103) 

(it coincides with Q-^-b (101)), and its domain is 

2^(aU) = {*V'.,B(r)}, (104) 



^Formally, a]^ ^ = a-^-B , therefore, again formally, 
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where ^,iJji,,B are absolutely continuous (inside R+) square integrable functions allowing 
the representation 



(105) 



1 di a_^^_B (0* 'P>y,B iO + • 
J ro 

Here, = ai^B*'^i',B, the image of on ^i/j^^b , is square integrable, ^(f)^^B £ 

L2 {R+) ; c^,B is a constant restricted by the requirement that *ipu,B G L2 (-R+) ; > 
is chosen appropriately depending on the value of u and the sign of i?. Of course, the 
same is true for al^_^ : we must simply make the substitutions u —u, B ^ —B in 
(103-105). 

At this point, we depart from the general procedure used in Sec. 3: then finding the 
deficiency subspaces and etc. Instead, we determine the "asymmetry" of evaluating 
the difference 

A = (,^', nt,^) - ( nt,^'„ m) (106) 

for any and**' belonging to Vi^Tl^). Using the form of G T>{Yl'^), see (102), (104), 
(105), 



and the form (102) of 11^^, we obtain (we here omit the subscripts * and B as irrelevant, 
which becomes clear below) 



/ dr ij'^ (r) al.ijl, (r) + ij% (r) ali^l (r) 



dr 



al,i)% (r) i^l (r) + al^)'^ (r) (r) 



Integrating by parts in the first integral, we find 

A = - 



{r)il}\ (r) + 11^% {r)il^l (r) 



i.e., A is determined by the asymptotic behavior of (107) at the boundaries, as r — > 00 

and r — > 0. The asymptotic behavior for ^V'l'.-B (?") and ^■.ip-u,-B (^) can be estimated using 
the representation (105) with the appropriate choice of Tq and estimating the integral 
term in (105) via the Cauchy-Bounjakowsky inequality. 

For example, we estimate the behavior of *ipu,B {r) at infinity, as r ^ 00, in the case 
S > 0. In this case, it is convenient to take ro = 00, such that 



*'^u,B (r) = i(7^,B (r 

The Cauchy-Bounjakowsky inequality yields 

eBe 



d^ (T^^^^B (0* (t>v,B (0 + Cz.,B 



(108) 



X 



der'^exp 



di\.<i>u,B (or 



< 



der'^exp 



eBe 



nl/2 



1/2 



-I/-1/2 

'eW 



exp 



eBr"^ 



1 + 



e(r) , (109) 
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where e (r) — as r ^ oo. It follows from (109) and the condition ^VVs ^ -^2 that 
c,^ B in (108) must be zero and, consequently, ^ipu,B {t^) vanishes as r — >• cxo, faster than 
r~^/^, independently from u. The same is true for B < 0, which is estabhshed via the 
Cauchy-Bounjakowsky inequality with ro < 00. 

The asymptotic estimates for (r) as r — > depend on the value of u. The result 

is: 

^ >^ ■ \*'4^u,B (r)| < c^r^/^ , 

^^t;- I *i^i/2,B (r) I < ci/2 r^/^ In — , 
z ro 

< < ^ : (r) = c,,B (rM)" + o (r^/') , 

independently from B; the factor M is introduced by the dimensional reasons. 

We conclude that for \u\ > 1/2, the two-spinors (r) (107) vanish both at infinity 
and at the origin, whereas for < 1/2 , the two-spinors (r) vanish at infinity, but are 
generally singular as r ^ : 

where c^'^ are arbitrary constants. With this estimates in hand, we find 

A- I 0,H>l/2 . . 

independently of S (we therefore omit the subscripts * and B as irrelevant) . 

We thus find that if > 1/2 , i.e., if Z 7^ or Z = and /i = 0, the operator is 
symmetric, 11''' C (n''^)\ and therefore self-adjoint, 11''' = (n^')^ it is sufficient to take 

the inverse inclusion (ll^)^ C II'I' into account, the standard relation for any symmetric 
operator (which is the consequence of the general relations A (1 B —>■ B^ Q for any 
densely defined operators). Thus, 11"'" is a self-adjoint extension of 11. 

This impUes that 11 is essentially self-adjoint, i.e., its closure 11 is self-adjoint, 11 = Il\ 
in addition 11 = 11^ and is a unique self-adjoint extension of 11. We recall the standard 
arguments. As is well known, H = (W) and H = Using the established self- 

adjointness of 11^, 11^ = (n^) , we obtain the chain of equalities 11 = (n^) = 11^ = 11 , 
which proves the first two assertions. There are no other self-adjoint extensions because 
any such extension H^"^* = (H^^^*)^ must satisfy the relation 11 C H^^* C W, H^^* must be 
the extension of 11 and the restriction of 11^^, but because 11 = H^^, there are "no place " 
for another self- adjoint extension. 

The situation is nontrivial in the case \u\ < 1/2, i.e. if Z = and 7^ 0, where the 
operator is nonsymmetric, which implies that 11 is only symmetric: 11 = (n^)^ C = 

is a strict inclusion. This relation allows finding 11. The inclusion 11 C 11^ implies that 
P (n) C P (nt), i. e. e X) (n) has the representations hke (105) and (110), 
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and vanish at infinity. Then the equahty 11 = (ll^") C 11''' implies, by the definition of 

{W) \ that ^' e P (IT) iff the difference A (106) vanishes for e P (IT) cV (W) and 
any ^ e P (n^)- According to (111), this gives 

+ c^c^ = Vc^ 

whence c'J = c',^^ = 0. H is thus defined as the restriction of to the domain V (ll) of 
the functions belonging to T> (H^), but vanishing at the origin. 

"We must seek the self-adjoint extensions H^^* = (11'^^*)^ of IT. These, as was said just 
before, "must he between" 11 and 11''', 11 C U^^^ = (ll"^^*)^ C 11'''. We start with finding 
the nontrivial (not coinciding with 11) symmetric extensions n*^™ of 11, 11 C n*'^™ C 
(jjsymy C and then show that these n*^™ solve the problem. The last inclusions 
allows repeating for n*^'" the previous arguments for n except that now the difference A 
(106) must vanish for any ^, ^' e P (n^?'"') C V (n^). According to (111) this gives the 
equation for cj,, ci^ in * (110) and c'J, c% in (112), 

^c%+^4 = 0, (113) 

Of course, if ^, ^' E V {U) C D(n^f™), i.e. vanish at the origin, Eq. (113) holds. 
To satisfy (113) it is sufficient that only ^ or ^' belong to V (IT). But D (n*^^™) must 
contain functions nonvanishing at the origin. Let be such a function and, for example, 
let cl^ ^ 0. For = Eq. (113) becomes 

+ ~^(^ = 

whence 

^ = -=^ = iA, AeM. 

Then for any ^' G D (11^^™) nonvanishing at the origin and this fixed ^f, Eq. (113) 
becomes 

^ + iA^=0, 

whence c'^^ ^ and dlld\ — iK, with the same A for all e I> (n*^"*) nonvanishing at 
the origin. 

The case c^^, = 0, 7^ is considered similarly. It is formally covered by the case 
A = ±cxo, where A = +00 and A = —00 are equivalent (both simply mean c^^ = 0, 
7^ 0). We thus obtain the asymptotic boundary conditions 

at the origin as r — > 0, that define a one-parameter family {n^j of all nontrivial sym- 
metric extensions of 11. 

It remains to show that these extensions 11 ^ are really self-adjoint, 11^ = (n^)^- For 

this, it is necessary to evaluate (11"^)^ Using the arguments similar to the previous ones, 

we conclude that ^' e P (^(n^)^j iff the difference A (106) vanish for ^' eV (^(H^)^) C 

V (nt) and any ^ e D (n^) C P (H'') , which yields 

^ + iA?2;=0. 
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If c'J and c'^^ are not equal to zero, this gives 

But this means that V (^(n^)^) C V (n^), i.e. (n^)^ C n^, the inclusion inverse to the 

initial C (n^)^ whence = (^^)^ which proves the final statement. 

We thus found the whole set of the self-adjoint extensions of the initial symmetric oper- 
ator n (100) with T> (n) =T>®T>. This is a one -parameter family {H'^}, —oo < A < cxo, 

whose each member 11^ = (hl'^)^ is defined by the asymptotic boundary condition (114) 
at the origin on * e P (n^j c V (nt). Because A = ±oo are equivalent, this family is 
homeomorphic to a circle U (1), not an open line R. The obtained self-adjoint asymptotic 
boundary conditions (114) evidently coincide with the self-adjoint asymptotic boundary 
conditions (45) (up to the common factor M^/^) obtained by the general method. As a 
by product, we find that the deficiency indices of 11 are (1, 1) in the case < 1/2. 

We conclude with the remark that those who is well aware of special functions (to 
effectively determine the deficiency subspaces) may prefer the general method (we note 
that, in fact, only the asymptotic behavior at the origin of the corresponding functions is 
need). But in any case, the evaluation of the adjoint operator, 11"'' or h\ is imperative in 
order to correctly determine the domain of the self-adjoint extensions. 

6 Concluding remarks 

We have studied the solutions of the Dirac equation with the magnetic-solenoid field in 
2-1-1 and 3-1-1 dimensions in detail. In the general case, there no simple relations between 
the solutions in 2 + 1 and 3 + 1 dimensions. However, wc have demonstrated that the solu- 
tions in 3 + 1 dimensions with special spin quantum numbers can be constructed directly 
based of the solutions in 2 + 1 dimensions. For this, we must choose the 2;-component 
of the polarization pseudovector as the spin operator in 3 -|- 1 dimensions. This is a 
new result not only for the magnetic-solenoid field background, but also for the pure AB 
field. The choice of as the spin operator is convenient from different standpoints. For 
example, the solutions with arbitrary momentum arc the eigenvectors of the operator 
S^. This allows explicitly separating the spin and coordinate variables in 3 -t- 1 dimen- 
sions and reducing the problem to the problem of the self-adjoint extension for the radial 
Hamiltonian only. Moreover, the boundary conditions in such a representation do not vio- 
late the translation invariancc along the natural direction which is the magnetic-solenoid 
field direction. Using von Neumann's theory of the self-adjoint extensions of symmetric 
operators, we have constructed the self-adjoint extensions of the Dirac Hamiltonian with 
the magnetic-solenoid field and obtained a one-parameter family and a two-parameter 
family of admissible self-adjoint boundary conditions in respective 2-1-1 dimensions and 
3 + 1 dimensions. The complete orthonormal sets of solutions thus have been found. 
We have determined the energy spectra dependent on the extension parameter for 
different self-adjoint extensions. In addition, we have for the first time described the 
solutions of the Dirac equation with the regularized magnetic-solenoid field in detail. We 
have considered an arbitrary magnetic field distribution inside the finite-radius solenoid 
and shown that the extension parameters © = — sgn($)7r/2 in 2 + 1 dimensions and 
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0_i_i = 0_i = — sgn($)7r/2 in 3 + 1 dimensions correspond to the limiting case R 
of the regularized magnetic-solenoid field. The finite radius solenoid regularization was 
also considered for the spinless particle case. It was demonstrated that in contrast to the 
spinning particle case, the corresponding (as it! — > 0) radial functions are regular for the 
arbitrary magnetic field inside the solenoid. 
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A Appendix 

l.The Laguerre function In,m{x) is defined by 



In,M - \ l\lT\v7^~'''^\ '^'"''"^^-'^^ n-m + l;x). (115) 
y r (1 -I- m) r (1 -I- n — m) 

Here, $ (a, b; x) is the confiuent hypergeometric function in a standard definition (see [38], 
9.210). Let m be a non- negative integer number; then the Laguerre function is related to 
the Laguerre polynomials L^{x) ([38], 8.970, 8.972.1) by 



Im+aA^) = y T{m + a + lf ^^'^"^'^'^^^^ ' ^^^^^ 
1 d"^ 

L° (x) = — -e^x-°^e-^x'"+" . (117) 



Using the well-known properties of the confiuent hypergeometric function ( [38], 9.212; 
9.213; 9.216), we can easily obtain the following relations for the Laguerre functions: 

Isjxin + l)4+i,„(j;) = (n - m x)ln,m{x) - 2x1'^ J^x) , (118) 

2v^(mTl)4,m+i(2;) ^{n-m- x)In,m{x) + 2xl'^^^{x) , (119) 

2^/xnIn~i,m{x) ^{n-m + x)In,m{x) + 2xl'^ ,^{x) , (120) 

2^/xmIn,m-i{x) ^{n-m- x)In,m{x) - 2xl'„,^{x) . (121) 

Using the properties of the confiuent hypergeometric function, we obtain the representa- 
tion 

, , r(l + n) exp(a;/2) ^-^ ^ , 

^"-<^) = V r(i + J)r(i+;-i.) -^^'^ + „, 1 + „ - (122) 

and the relation ([38], 9.214) 

In,m{^) = i-ir-'^ImA^), u-m integer . (123) 
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The functions Ia+m,m{^) satisfy the orthonormaUty condition, 



^0 



which follows from the corresponding properties of the Laguerre polynomials ( 
7.414.3). The set of the Laguerre functions 

Ia+m,m{x)-, m = 0, 1, 2... , a > -1 

is complete in the space of square integrable functions on the half-line {x > 0), 

oo 

Ia+m,m{x)Ia+m,m{y) ^ S {x - y) . (125) 

m=0 

2. The function V'A.al^) is even with respect to the index a, 

1pX,a (x) = 1pX-a (x) ■ (126) 

It can be expressed via the confluent hypergeometric functions, 



^ ^ X,l + a;x 



r(i^-A) V 2 



* — X,l- a;x 



(127) 



'r(i±^-A) 

or, using (115), via the Laguerre functions, 

/ X \/r (1 + //) r (1 + m) , , , , 

il>Xa (x) = — ; ^ ismrmlnm (x) — sin m7r/^ „ (x)) , 

^ ' ^ ^ sin (n - m) TT ^ ' ^ ^ rn,ny JJ , 

1 — a 1 + CK 
a — n — m, 2X — 1 + n + m, n — X — , m — X — . (128) 

The relations 

/\ ^ : /\ l + a — 2A, 

^\,a [X) = Vxi/jx-^a-l {X) + 2 ' 

1pX,a (x) = Vxipx-la+l (^) + ' ^-^V'A-l,a (x) , 

2a;^l,a (x) = (2A - 1 - x) ^A,a (x) + i (2A - 1 - a) (2A - 1 + a) V^A-i,a (a;) , 
22;'0l,a (2;) = (a - x) V'A.a (a^) + (2A - 1 - a) Vxipx-^a+i (x) 

= (x-2A-l)V'A,a-2V'A+l,a (129) 

are valid for the functions ip\^a (x) . The direct consequence of these relations is 

Aa1pX,a (x) = ^^-^-^V'A-i,a-l i^) > ^a'^X-^a-l i^) = V'A.a (x) , 
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Using the well-known asymptotic behavior of the Whittaker function ([38], 9.227), we 
have 



V'A,a (x) r^X^ 26 2 , X ^ OO; VA,a (x) ~ 



r(l«l) 



_|a| 



x-^ , a 7^0, x~0. (131) 



The function ip\^a (x) is correctly defined and infinitely differentiable for < x < oo and 
for any complex A, a. In this respect , we note that the Laguerre function is not defined 
for negative integer n and m. In the particular cases where one of the numbers n or m 
is non-negative and integer, the function ipx^a (x) coincides (up to a constant factor) with 
the Laguerre function. 

According to (131), the functions ipx^a (x) are square integrable on the interval < a; < 
oo whenever \a\ < 1. This is not true for \a\ > 1. The corresponding integrals for a 7^ 
are calculated (see 7.611 in ([38] )): 

^A,a {x) %lJx',a {x) dx 



l-a-2X' 



(A' — A) sin an 
l + a-2X 



2 

|q;| < 1 , 



1 - « - 2A 



00 

/ 



\tp\,a (x) I'^dx 



71 ^(^±^)-^(^^) 

sinaTT r (i±^) r (i^^f^) ' 



lal < 1 . 



(132) 
(133) 



Here, '^(x) is the logarithmic derivative of the F-function ( [38], 8.360). In the general 
case, the functions ipxa (x) and ipy a{x), A' 7^ A, are not orthogonal, as it follows from 
(132). 



B Appendix 

We here present modifications of some foregoing formulas for the case B < 0. 

1. The spectrum uj corresponding to the functions 4>m,i,a{'r) is 

r 27(m-/ + l-/i), /-(l + a)/2<0 
"^"1 27(m+(l-a)/2), Z-(l + a)/2>0 ' 

and the spectrum uj corresponding to the functions 

^^['^1(^ + ^-1^)^ '^^-l . (135) 

[ 27m, (7 = 1 ^ ' 

These expressions are the modifications of Eqs. (18) and (19) for the case S < 0. 

2. We consider the spinors ui (r) satisfying (20). In the case a; = 0, these are 

= ( 1""'-''^ ) . i > 1; = ( f'*--' ) , i = 0. (136) 
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In the case 7^ 0, these are 



«i:.:,.W=(f;^»fj,,,).' = 0.^ = 27(m+l), 

= ( t7«'!'(r) ) . ' = . = 2t + 1 - /<) . (137) 

The presented expressions are the modifications of Eqs. (21) and (22) for the case B < 0. 

3. In the case a; = 0, the only positive energy solutions (particles) of Eq. (10) are 
possible. These solutions coincide with the corresponding spinors m up to a normalization 
constant: 

+Mr) = NuoAr), 1>1; +i^i\r) = Nul\r), 1^0. (138) 

Thus, the particles have the rest energy level, whereas the energy spectrum of antiparticles 
begins from _e = —^/M"^ + 27. 

4. For the case S < 0, the relations for the irregular spinors u^^^c{J) similar to the 
relations (27) are 



Ilii^ _i(r) = i^f2riu^;y{T), ^u^^t) = -i—=u^_x{r) ■ (139) 
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